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ABSTRACT 

We report the results of the analysis of the extended soft emission discovered in the Chandra ACIS 
pointing at the merging system NGC4038/39 (the Antennae). We present a 'multi-color' X-ray image 
that suggests both extensive absorption by the dust in this system, peaking in the contact region, as 
well as variations in the temperature of different emitting regions of the hot interstellar medium (ISM). 
Spectral fits to multi-component thermal emission models confirm this picture and give a first evaluation 
of the parameters of the hot plasma. We compare the diffuse X-ray emission with radio continuum 
(6cm), HI, CO, and Ha images to take a first look at the multi-phase ISM of the Antennae galaxies. 
We find that the hot (X-ray) and cold (CO) gas have comparable thermal pressures in the two nuclear 
regions. We also conclude that the displacement between the peak of the diffuse X-ray emission in the 
north of the galaxy system, towards the inner regions of the northern spiral arm (as defined by Ha, radio 
continuum and HI), could result from ram pressure of infalling HI clouds. 

Subject headings: galaxies: peculiar — galaxies: individual (NGC 4038/39, Antennae) — galaxies: 
interactions — X-rays: galaxies 



1. INTRODUCTION 

This is the fifth in a series of papers based on the first 
Chandra ACIS-S (Garmire 1997; Weisskopf et al 2000) ob- 
servation of the Antennae galaxies ( NGC4038/39). The 
first paper in this series gives a general picture of the over- 
all results of this 72.2 ks. observation (ObsID 315: Fab- 
biano et al 2001 - Paper I). The following three papers 
(Zezas et al 2002a, b - Paper II, HI; Zezas & Fabbiano 
2002 - Paper IV) discuss the properties of the point-like 
source population detected in this merging system. Here, 
we concentrate on the analysis of the complex diffuse soft 
emission reported in Paper I and on its importance for our 
understanding of the hot ISM of these galaxies. 

The Antennae (dynamically modeled by Toomre & 
Toomre 1972 and Barnes 1988) are the nearest pair of col- 
liding galaxies involved in a major merger (D = 19 Mpc for 
Ho = 75 km s^^ Mpc~^ ^). Hence, this system provides a 
unique opportunity for getting the most detailed look pos- 
sible at the consequences of a galaxy merger as evidenced 
by induced star formation and the conditions in the ISM. 
Each of the two colliding disks shows rings of giant H II 
regions and bright stellar knots with luminosities up to 
My — 16 (Rubin et al 1970), which are resolved with 
the Hubble Space Telescope into typically about a dozen 
young star clusters (Whitmore & Schweizer 1995, Whit- 
more et al 1999). These knots coincide with the peaks of 
Ha, 2.2/x, and 6-cm radio-continuum emission (Amram et 
al 1992; Stanford et al 1990; Neff & Ulvestad 2000), indi- 
cating an intensity of star formation in each single region 
exceeding that observed in 30 Doradus. CO aperture syn- 



thesis maps reveal major concentrations of molecular gas, 
including ~2.7x10^Mq in the region where the two disks 
overlap (Stanford et al 1990; Wilson et al 2000). A recent 
K-band study derives ages for the star clusters ranging 
from 4 to 13 Myr, and measures high values of extinction 
{Av ~ 0.7 - 4.3 mag; Mengel et al 2000). 

The presence of an abundant hot ISM in the Anten- 
nae was originally suggested by the first Einstein observa- 
tions of this system (Fabbiano & Trinchieri 1983), and 
has been revisited with all following major X-ray tele- 
scopes (ROSAT: Read et al 1995, Fabbiano et al 1997; 
and ASCA: Sansom et al 1996). The ROSAT HRI im- 
age was used in a recent multiwavelength study of the 
Antennae by Zhang et al (2001), which suggests feedback 
between star clusters and the interstellar medium (ISM). 
The Chandra observation gives us the means to study this 
hot ISM in unprecedented detail. The spatial resolution 
of Chandra is at least 10 times superior to that of any 
other X-ray observatory, allowing us to resolve the emis- 
sion on physical sizes of ~ 47 pc for D=19 Mpc. With 
this resolution we have detected individual X-ray sources 
down to Lx ~ lO^^ergs (Paper II) and we can image 
in detail the spatial properties of more extended emission 
regions (see Paper I). At the same time, ACIS allows us 
to study for the first time the X-ray spectral properties 
of these emission regions, providing additional important 
constraints on their nature. 

Our initial look at these Chandra data (Paper I) revealed 
a complex, diffuse and soft emission component responsi- 
ble for about half of the detected X-ray photons from the 



^ We have adopted this distance for conformity with most recent work on the Antennae; in Papers I-III, a distance of 29 Mpc, 
50 km s~^ Mpc~^ was used; in Paper IV we discuss our results using both D = 19 and 29 Mpc 
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two merging galaxies. The X-ray spectrum of this com- 
ponent exhibits emission Hnes, pointing to optically thin 
gaseous emission. In the present paper, we explore fur- 
ther the spectral/spatial properties of this emission and 
its relation to the multi-phase ISM: in §2 we present a 
'true-color' X-ray image of the extended emission (i.e., an 
image where the color scale represents the photon energy), 
in §3 we present a spectral fitting analysis of these data, 
and in §4 we describe new 4 and 6cm radio images derived 
from archival radio data, matched to the effective resolu- 
tion of the diffuse X-ray emission. In §5 we discuss the 
physical properties of the hot ISM and we compare the 
diffuse X-ray emission with Ha, radio continuum (6cm), 
HI and CO data. 

2. THE TRUE-COLOR X-RAY IMAGE 

Chandra was pointed to NGC4038/39 on December 1, 
1999, for 72.2 ks with the back-illuminated ACIS-S3 CCD 

chip at the focus. The observations and data reduction 
and screening are described in Paper I. We used CIAO 
2.2.1 and other software developed for the purpose in the 
present data analysis. 

We have used these data to produce a true-color X-ray 
image of the Antennae that would be sensitive to spectral 
differences in the soft diffuse emission. The true-color X- 
ray image of Paper I showed that the diffuse emission is 
soft (mostly at energies below 2 keV) and thus differenti- 
ated it from the harder broad-band emission of the point- 
like sources, but did not have the resolution to show spec- 
tral differences in this soft emission. To obtain this resolu- 
tion, we cut the data into energy bands of 0.3-0.65, 0.65- 
1.5, and 1.5-6 kcV. These bands were chosen using thermal 
plasma models as a reference: the low energy band does 
not include the Fe-L blend seen around 1 keV, whereas the 
high energy band should include only the hottest (or most 
highly absorbed) material. The resulting images were then 
adaptivcly smoothed using the tool csm,ooth with a min- 
imal significance (S/N) of 3.5 and maximum significance 
(S/N) of 5.0. This tool uses different smoothing scales 
to enhance low-surface-brightness extended features, while 
retaining the resolution allowed by the areas of the image 
with the highest signal-to-noise ratios. 

Monochromatic exposure maps were created for each of 
the bands at 0.5, 1, and 3 keV, and smoothed using the 
same smoothing scales as the images. The images were 
then divided by the corresponding exposure map and com- 
bined using dmimg2jpg to create a "true-color" image, with 
a red to blue progression to represent increasing photon 
energy (Fig. 1). 

To derive a true-color image of the diffuse emission, we 
then proceeded to identify and remove point-like sources. 
A source list was first generated using the wavdetect tool, 
using scales of 1, 1.4, 2, 2.8 and 4 pixels; the initial list was 
then culled to include only obvious point sources. These 
arc sources 1, 2, 3, 4, 11, 12, 13, 14, 16, 18, 19, 20, 21, 23, 
29, 31, 32, 33, 34, 35, 37, 38, 41, 42, 43, 44, 46, 47, 49, 
from the list of Paper II. Background regions were created 
to best represent the local diffuse emission. The tool dm- 
filth (dmFillTheHole) was run in each of the three spectral 
band images to fill in the source regions with pixel values 
sampled from the Poisson distribution whose mean was 
that of pixel values in the corresponding background re- 



gion. The same procedure described above of smoothing, 
exposure correction and combining of the three resulting 
images into a true-color image was then followed. The 
result is shown in fig. 2. 

While figs. 1, 2 do not give a quantitative assessment of 
the spectral properties, because the color scale is not linear 
and has been enhanced to show faint features, it is clear 
from them that there are spectral variations throughout. 
In particular, the area of impact between the two galaxies 
[centered approximately on RA(J2000)=12''01™55* and 
Dec(J2000)=— 18°52 50 ] releases only high energy pho- 
tons, as it would be expected in the presence of large 
amounts of obscuring material. The nuclear regions, as 
already suggested in Paper I, include harder emission 
('white' in figs 1, 2) than other relatively unobscured re- 
gions of the ISM. 

3. SPECTRAL ANALYSIS OF THE DIFFUSE X-RAY 

EMISSION 

To derive quantitative information on the spectral differ- 
ences suggested by the color maps, we performed a spectral 
analysis of the diffuse emission. We screened the data for 
all the point-like sources reported in Paper II. We retained 
the sources with a strong extended component. These are 
sources 5, 6, 7, 10, 25 of Paper II. Fig. 3a shows a raw 
image of the Antennae together with the excised source 
regions. We did not correct for the amount of diffuse emis- 
sion which was also excised together with the excluded 
sources. This does not affect the parameters of the fit 
(temperature, Nh, photon index, abundance) but results 
in lower normalizations (or emission measure in the case 
of a thermal model) . 

To derive quantitative information on the varying spec- 
tral properties suggested by the color image, we first di- 
vided the area in a grid, so to perform a systematic spec- 
tral analysis from which we could derive a first-cut spectral 
map of the hot ISM. In order to have at least 100 in each 
extracted spectrum (before background removal), we used 
regions 50 x 50 pixels (24.5" x 24.5") in size. Using fig. 2 
as a guide, we then focused on the regions with the largest 
number of counts (highest signal to noise) to perform a 
more in depth analysis. Both grid and individual regions 
(ellipses) are shown superimposed on the color map of the 
diffuse emission in fig. 3b. 

For each region (either square in the grid, or ellipse) 
we extracted photon invariant (PI) spectra and created 
weighted response matrices (RMFs) and ancillary response 
files (ARFs) using the acisspec script of CIAO v2.1. This 
script invokes the dmextract tool to extract the spectrum 
and then the mkwrmf and mkwarf tools in order to cre- 
ate area weighted RMF and ARF files for each region, 
since they span more than one Fits Embedded Function 
(FEF) regions. The spectra were then grouped to have at 
least 15 counts per bin in order to use statistic for the 
fits. From the spectrum of the diffuse emission presented 
in Paper I, is clear that most photons are below 5 keV. 
Therefore in order to minimize contamination by the back- 
ground which dominates above 7 keV we ignored any data 
outside the 0.3-7.0 keV energy range. For the spectral fits 
we used the XSPEC vlO.O spectral fitting package (Arnaud 
et al. 1996). All the cited errors for the spectral param- 
eters are for the 90% confidence level for one interesting 
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parameter, unless otherwise stated. The background was 
obtained from source free regions outside the galaxy. 

3.1. Spectral 'Map' 

First we attempted fits with an absorbed single ther- 
mal plasma model (RS model; Raymond & Smith 1977) 
and solar abundances (the abundance table by Anders 
& Grevesse 1989 was used). For these fits the absorb- 
ing column density was free to vary above the Galactic 
= 3.4 X 10^0 cm-2 (Stark ct al. 1992). Typical best- 
fit kT is betwcccn 0.1 kcV and 0.9 kcV and best-fit Nr 
can reach as high as ~ 8 x 10^^ cm~^ . However, these fits 
were unacceptable in all but 6 cases (reduced greater 
than 1.4). Letting the metal abundance vary as a free pa- 
rameter in the fits (with the relative abundances fixed at 
the solar ratio) increased the number of fits with accept- 
able x^, but the fit is still unacceptable for the majority 
of the spectra. The values of the abundances suggested by 
the fits were uniformly subsolar (^ 0.1). 

Since a single thermal component model is clearly not an 
adequate representation of the data, we introduced a sec- 
ond thermal component with a different absorbing column 
density [RS-(RS) model; here and in the following mod- 
els the parentheses indicate the absorbed component] . We 
first again assumed solar abundances, and we then let the 
abundance vary in both components. In both cases, the 
overall N/f was free to vary above the Galactic value. 

The two-component, solar abundance, model provided 
a statistically significant improvement in the fit (above 
the 99.9% confidence level based on an F-test for two ad- 
ditional parameters; Bevington & Robinson, 1992) in all 
but 6 cases, over the single component, solar abundance 
model. In most cases the overall absorbing column is very 
close to the Galactic Nh and the low temperature com- 
ponent has a typical kT of 0.2-0.3 keV. The second com- 
ponent is typically seen through a larger column density 
( Nh^ 2 — 8 X 10^^ cm~^ ), and has a higher temperature 
than the unabsorbed component (kT~ 0.7 — 0.8 keV). The 
results from these fits are presented in Table 1: Golumn 
(1) gives the data extraction region following the nota- 
tion of fig. 3b, Column (2) gives the net number of counts 
in each region (together with the statistical error). Col- 
umn (3) gives the overall Nh in units of 10^^ cm~^ 
, Column (4) gives the temperature for the first compo- 
nent, while Columns (5) and (6) give the column density 
and the temperature for the second (absorbed) compo- 
nent. Column (7) gives the ratio of the emission measures 
of the two best-fit thermal components. Column (8) gives 
the and the number of degrees of freedom (dof). For 
the regions where a single component model gives an ac- 
ceptable fit only the parameters for this model are listed. 
Fig. 4 presents the 99% confidence contours for the two 
interesting parameters Nh and kT. Red and blue con- 
tours identify the low and high temperature components 
respectively. The crosses mark the best fit values for each 
component; if the crosses are not shown the best fit values 
fall outside the range of the plot. Here, as in all the spec- 
tral fit tables, the confidence regions are calculated usind 
the Ax^ method, therefore they are a strong constraint 
only in the case of acceptable fits (i.e. when the model 
is a reliable representation of the data). Regions 2C, 2E, 
3E and 5G have unacceptable x^ in the solar abundance 



fit. These areas arc part of the luminous emission regions 
addressed in greater detail in §3.2. In two instances in 
Table 1 (and Table 2) the best-fit kT of the low temper- 
ature component is below 0.1 keV (i.e. outside the range 
of ACIS). However, in these cases the upper bound is typ- 
ically well within the instrumental range. These results 
should be considered conservatively as upper limits on the 
temperature. 

The results of the fits with the abundance as a free pa- 
rameter in each thermal component are listed in Table 2. 
Although the x^ are slightly improved, the fits are still 
not acceptable in the bins with a larger number of counts. 
The best-fit temperatures, although not identical, are in 
the general range of those obtained for the solar abun- 
dance fits. The abundance is either ill-defined or subsolar 
[abundances are listed in Cols. (4) and (7)]. Whenever 
a minimum is found, we list it with the errors. In some 
cases, a minimum is foimd, but the confidence region ex- 
tends outside one of the fit boundaries (0-5 for Z; 0-64 keV 
for kT; > 3.4 x 10^° cm~^, the line of sight value, for the 
total Nif). In these cases, we list both the value of the pa- 
rameter in question at the minimum x^ and the bound that 
can be determined from the explored grid. If no minimum 
is found, we give the bound calculated on the basis of the 
fitting grid boundary where the lowest value of x^ is regis- 
tered. In the case of the abundance, large values typically 
occur in unconstrained cases, either because a minimum 
was not found in the x^ distribution, or, if a minimum 
was found, because the error is of the order of magnitude 
of the best fit value. Some values of Njy returned by the 
fit (e.g. for regions 5-D and 5-E) are extremely large, but 
also have comparably larger errors, indicating that they 
are likely to be a consequence of poor statistics. In the 
case of region 5-E, the extreme value of Nh is not found 
when other models are adopted (e.g. see Table 1). Simi- 
larly, the ratios of best-fit emission measures (both Table 1 
and 2) show a few extreme values. This occurs for regions 
2-A, 5-E (both tables), and 2-E, 5-D (Table 2), indicating 
that one of the two components is likely to be an artifact 
of the fit. 

We also fitted the data with a composite thermal plasma 
-|- power-law model (RS-I-PO). The rationale for this model 
is that there could be a significant unresolved X-ray binary 
component to the emission, with fluxes below the detec- 
tion threshold of our data (see Paper II; the corresponding 
luminosity is near lO'^^ergs s^"'^). This model gave consis- 
tently worse x^ than the two-thermal-plasma model. 

In summary, we conclude that this first fitting exercise 
strongly suggests a complex temperature distribution in 
the ISM, at least in all the regions with higher count statis- 
tics. 

3.2. Four Luminous Emission Regions 

Although the spectral grid gives a quantified version 
of the color image, showing spectral variability through- 
out the soft emission, and suggesting a multi-temperature 
ISM, the division of the data in a grid cuts some of the 
features of the ISM in a way that it is not optimal. We 
have therefore selected four regions, identified in fig. 3b, for 
further study. These regions are associated with the two 
nuclei and the chain of hotspots to the West of the north- 
ern nucleus (NGC 4038). In fig. 3b they are marked as N 
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for the northern nucleus (NGC 4038), S for the southern 
nucleus (NGC 4039), and Rl and R2 for the two hotspots 
to the West of the northern nucleus. The emission peak 
in the S region was excised, because it corresponds to a 
point-like source (Paper II). 

Wc repeated the fitting procedures used for the binned 
array, using again on these regions the RS+(RS) and 
RS+(PO) model, with both solar abundance and with the 
abundance as a free parameter. In the assumption of so- 
lar abundance, neither model gave an acceptable fit to the 
spectrum of any region but R2, which is also the region 
with the smallest number of spectral counts (see Table 3). 
We therefore did not pursue further analysis for R2. The 
results from these fits are summarized in Tables 3 and 4. 

Again, the introduction of the abundance as a free pa- 
rameter leads to improved, but still not good, fits. Pre- 
ferred values of the abundance are subsolar. This results is 
understandable if wc look at fig. 5, that shows the best fit 
RS-I-(RS) solar abundance model together with the residu- 
als for these four regions. The model soft component tends 
to have a more pronounced peak in the region below 1 keV 
than the data. This is the regions of more prominent Fe-L 
emission. With the exclusion of the RS+[PO] model for 
the two nuclei, the improvement of the fit over the solar 
abundance models was statistically significant above the 
99% confidence level based on an F-test. 

Fits to 3-component models, with fixed solar abundance 
(Table 5), give improved at a level similar to those of 
the free abundance 2-component fits. These models con- 
sist of two thermal RS components and either a power-law 
(PO) or a third RS component. One of the RS components 
is either absorbed by the same column as the power-law 
[RS-|-(RS+PO) model] or just by the foreground column 
[RS-I-RS-I-(PO) model]. Similarly, in the 3-RS compo- 
nent models, we also explored the possibility that the sec- 
ond thermal component is associated either with the fore- 
ground or with the obscured thermal component: these are 
noted as RS-hRS-h(RS) and RS-h(RS-hRS) respectively. 

All models gave similar x^, with the exception of the 
southern nucleus, where the 2-RS-l-PO models are slightly 
favored. The best fit temperatures for the two thermal 
components are kT~ 0.2 - 0.3 keV and kT~ 0.7 - 0.9 keV 
respectively, and the slope of the power-law is F ~ 3. 
These power-laws are steeper than expected from unre- 
solved X-ray binaries (where F ~ 2 are more the norm; 
e.g. Schulz 1999; Yokogawa ct al 2000). The results from 
the three thermal component models are also presented in 
Table 5. While two of the components are in the range 
found by the previous fits, the third component tends to 
be rather hard, and could be suggestive of unresolved hard 
sources. 

The range of models used are still not adequate to fit 
the spectra with the largest number of counts (southern 
nucleus and region Rl , in that order) . Wc further explored 
the worst-fit, highest counts southern nucleus, by allowing 
for the abundance to vary in the RS+(RS+PO) model. 
However, this did not help the fit (x^=95.8/63), and it 
led to totally unconstrained results for the power-law and 
the abundance of the soft component. We also divided the 



southern nucleus region in two sub-regions, one including 
the E-W bar-like feature (see fig. 2), but we did not get 
appreciably difi^erent results within statistics; the results 
are consistent with those of Table 5. 

These result suggest complex emission regions. The dif- 
ferent fits, both with solar abundance and with the abun- 
dance as a free parameter, suggest two principal tempera- 
tures of the hot ISM: ~ 0.2 - 0.3 keV and ~ 0.6 - 0.8 keV. 
The abundance is either ill-defined or subsolar. Subso- 
lar abundances were also reported by Sansom et al (1996) 
from the spectral analysis of the integrated ASCA spec- 
trum of The Antennae. At this point, we believe that 
one should be cautious in over-interpreting these possi- 
ble subsolar abundances, given that the complexity of the 
emission regions, and the relative low spectral resolution 
of the CCD data could conspire in giving erroneously lower 
abundances (see e.g. demonstration in Kim & Fabbiano 
2003: Weaver et al 2000; Matsushita et al 1997). More- 
over, the elemental abundances were fixed at the solar ra- 
tio, and this may not hold in a starburst regions with a 
prevalently young and evolving massive stellar population, 
where Iron may be under- abundant, because of the preva- 
lence of SN II. Finally, uncertainties in the line emission 
models themselves may conspire in giving an erroneously 
higher continuum, if the data are not of sufficient spectral 
resolution. Given the obvious spectral complexity of this 
ISM, and the relatively limited spatial resolution, resulting 
from the need to integrate the spectra over sizeable areas, 
we believe that delving deeper in the abundance problem 
is not warranted by these data. 

In Table 6 we give the emission measures for the fits of 
Table 4 (EM^) and power-law (PO) normalization for the 
three component models, together with the absorption- 
corrected luminosity in the 0.1 - 10.0 keV band for each 
component. 

4. LOW-RESOLUTION 4CM AND 6CM RADIO MAPS 

Since we do not sec sharp features in the hot ISM, pos- 
sibly because of the relatively low statistics of the X-ray 
data, we have re-analyzed radio wavelength observations 
to produce images of the diffuse radio emission, on scales 
similar to those of our true-color X-ray map (fig. 2) . 

VLA^ archival data (20, and 6cm) were used together 
with 6cm observations from Neff & Ulvestad (2000) to pro- 
duce images sensitive to a range of spatial scales. The stan- 
dard source 3C286 was used to set the flux density scale, 
and the source 1159-2148 was used as a phase calibrator. 
Each data set was edited, calibrated, imaged, and then 
self-calibrated using standard image processing techniques 
and the Astronomical Image Processing System (AIPS). 
Data sets observed using B1950 coordinates were precessed 
to the J2000 system. Data sets from different arrays were 
then combined and self-calibrated. Data weighting was ad- 
justed to produce images of the desired angular resolution. 
Image deconvolution (CLEANing) included the removal of 
four background radio sources within the telescope's pri- 
mary beam response. Because our final images are made 
from multiple data sets, some with unknown information 
about weather during the observations, we conservatively 



^ EM = J" n^dV where n is the density, V is the volume of the gas 

^ The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory (NRAO), which is a facility of the National Science 



Foundation operated under cooperative agreement by Associated Universities, Inc. 
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assume a flux calibration error of 5%. Absolute position 
errors of the brighter radio sources are less than 0.15". 20 
and 6cm images are shown in fig. 6, with a restoring beam 
size similar to that of the X-ray color map shown in fig. 2. 
The resulting images were low-pass filtered to remove com- 
pact radio sources (supernova remnants and HII regions), 
using a box 7.6 x 7.6 arcsec. The smoothed images were 
combined to form a spectral index image (spectral index 
a, where Si, oc using regions with S/N of 10 or better 
in the smoothed images (fig. 7). 

5. DISCUSSION 

5.1. Properties of the Hot ISM 

Fabbiano et al (2001) reported that the Chandra image 
of NGC4038/39 shows extended, soft emission features 
with dimensions ranging from a few arcseconds (few hun- 
dred parsecs) to the size of the entire galaxian system or 
larger. In §2 and §3 we have taken a closer look at this 
emission and at its spectral properties. Here we discuss 
the implications of these results for the properties of the 
hot ISM. 

We used the results of the spectral analysis to estimate 
the physical parameters of the ISM in the two nuclear re- 
gions and in the most intense well-defined emission region 
in the N-W arm (Rl, see fig 3b and Table 5), and to infer 
the supernova rate in these regions. We used the emis- 
sion measures from the fit to the RS + (RS + PO) model 
(see Table 6). In this model we assume a cooler external 
unabsorbed ISM layer surrounding hotter ISM and harder 
sources of emission embedded in a cooler /dusty absorbing 
medium. This picture agrees with the true-color X-ray 
map, and suggests a hotter inner superbubble, heated by 
early-type stellar winds and supernovae, with an outfiow- 
ing and cooling upper layer of ISM. 

Our estimates are given in Table 7, for both low-kT 
(subscript 1), and high-kT (subscript 2) thermal compo- 
nents (from EMi and EM2 of Table 6). For these cal- 
culations, we have assumed cylindrical emitting volumes 
of base given by the count extraction areas of Table 2. 
Although it is likely that hot absorbed emission originates 
from inner regions in a superbubble, with our data it is not 
possible to make an informed discrimination between the 
different geometries (even in projection on the sky) of the 
emitting regions. We assumed a height for the emitting 
cylinders of cither 200 pc or 1 kpc. The thinner cylinder 
would correspond to an emission region of the approximate 
depth of a spiral disk, while the thicker cylinder assumes a 
scale- height more consistent with an outflowing hot ISM. 
We give both estimates in Table 7. We estimated thermal 
energy and cooling times following Tucker (1975) and us- 
ing the intrinsic luminosities from Table 6. Thermal pres- 
sure in the assumption of electron and proton components 
in thermal equilibrium is p = 2nkT. 

Supernova rates were derived as in Heckman et al (1996; 
see also Fabbiano et al 1997 for a similar approach applied 
to the ROSAT HRI image of the Antennae), which as- 
sumes that most of the mechanical energy supplied to the 
ambient medium in the star-forming regions is from super- 
novae. Following Heckman et al (1996), we adopt the rela- 
tion for a superbubble: Eth = 1.4x 10^'^Lmec/t,43i7ergs s~^, 
whre Lmech,43 is the mechanical energy supplied to the 
ambient medium in the star forming region from super- 



novae, and tj is the age of the starburst in units of 10*" yr. 
Lmech = RsN X EsN^ whcrc RsN is the supernova rate 
and EsN = 10^^ ergs is the energy released in a supernova 
explosion. Prom this, we obtain an estimate of the super- 
nova rate Rsn, equating Eth to our estimates of Table 7. 

As shown by the X-ray/Ha comparison in Paper I, in 
the nuclear regions the hot and warm ISM distributions 
closely follow each other, suggesting that the two gases 
are intermingled, as for example in the gaseous outflows 
of the nearby starbiu'st nucleus AI82 (Watson et al 1984). 
In this case the X-ray emitting volume may well be a frac- 
tion {rj) of the one adopted here [t] = 1 was assumed 
in Table. 7). This would result in shorter cooling times 
(tc ~ T}^^^) and similarly smaller thermal energies and su- 
pernova rates, while the thermal pressure would increase 
as p ~ J7~^/^. 

The other region in Table 7 (Rl) was identified in Pa- 
per I as a superbubble, with a hot core surrounded by Ha 
filaments. As suggested in Paper I, Rl may be a more spec- 
tacular example of the superbubbles seen in more nearby 
galaxies such as the LMC and MlOl (e.g. Wang & Helfand 
1991; Williams & Chu 1995; see reviews in Chu 2000 and 
MacLow 2000). 

We derive cooling times in the 10^ — 10* yr range and 
total masses of hot ISM in the 10^ - 10^ for each 
region. For the two nuclear regions, where a comparison 
is possible, these hot gas masses are of order 1/1000 of 
the mass in cold molecular gas (Wilson et al 2000). The 
thermal pressures of the hot gas are of the order of a few 
X 10~^^ dyn cm~^. For the two nuclear regions, these pres- 
sures are remarkably similar to those that can be derived 
for the cold molecular clouds from the CO measurements 
of Zhu, Seaquist & Kuno (2003), suggesting a rough pres- 
sure equilibrium of the different phases of the ISM. Using 
the values of densities and temperatures from Zhu et al, 
we obtain cold gas pressures of 4.2x10^^^ dyn cm^^ and 
3.1x10"^^ dyn cm~^ for the northern and southern nu- 
clear regions, respectively. In both nuclei, for our assump- 
tion of gas volumes, the hotter X-ray components may 
have larger pressures, and create bubbles or even escape 
from the enveloping molecular clouds. However, we must 
remember that these estimates (both CO and X-rays) are 
only average values, depending on a number of assump- 
tions, including the emitting volume and the small scale 
properties of the clouds (density and temperature struc- 
tures) that cannot be measured at present. 

In these regions our estimates of the supernova rates 
are within factors of a few of those estimated by Neff & 
Ulvcstad (2000) from compact radio sources, on the basis 
of radio continuum observations at 6cm. We stress that 
the estimates in Table 7 are only indicative, given the un- 
certainties in the geometry of the emitting regions, and in 
the spectral models. If the 3-temperature RS-I-(RS-I-RS) 
model is used (see Table 5), the total emission measure 
would be larger for the northern nucleus. But even so, the 
estimates of the ISM parameters would not change much: 
this would result in a factor of < 2 increase in the mass of 
hot ISM. and correspondingly shorter cooling time. The 
estimates for the; southern nucleus would not change ap- 
preciably; in this case the fit returns two low-temperature 
components with temperature and EM similar to those 
from the RS-|-(RS-|-PO) model, and a higher temperature 
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component, that may be related to harder emission from 
supernova remnants and X-ray binaries. In the case of Rl, 
the results of this fit are probably unrealistic, requiring an 
extremely luminous, very low energy component. 

5.2. Large-Scale Morphology of the Multi-Wavelength 
Diffuse Emission 

Fig. 8 shows a comparison of the X-ray image of the dif- 
fuse emission (likely to be gas shock-heated by supcrnovac 
and stellar winds), with the HST WFPC-2 archival Ha 
image (Whitmorc ct al 1999), that traces photoionization 
of the ISM by luminous O stars, and with the low-pass fil- 
tered 6cm radio continuum VLA image of fig. 6, to which 
both the low-energy spectrum of the warm and hot ISM, 
as well as non-thermal contributions from embedded su- 
pernova remnants, contribute. Fig. 9 compares the diffuse 
X-ray emission of fig. 2 (here plotted in red) with CO (Wil- 
son et al 2000) and HI (Hibbard et al 2001) images. 

These comparisons give a fairly complete picture of the 
large scale distribution of the different components of the 
ISM, and also of the effect of the extinction on the emis- 
sion. 

The regions where young stars are being formed in dusty 
molecular clouds, as traced by the CO, are also promi- 
nent in the radio continuum, that is not affected by cold- 
ISM/dust absorption. The radio continuum appears to fol- 
low the Ha reasonably well, with the exception of the peak 
at the NE of the southern nucleus (the dusty star-forming 
'impact region', Mirabel et al 1998; Wilson et al 2000; 
Zhang et al 2001). The radio spectral index map (fig. 7) 
indicates three regions where thermal emission (a > 0.0)is 
present, probably from HII regions: the thickest part of 
the overlap region to the NE of NGC 4039 nucleus, another 
HII region complex further N and E of that, and a super- 
bubble structure at the southern edge of the arc extending 
SW from the NGC 4038 nucleus. Filamentary structures 
in the non-thermal emission are evident running roughly 
EW across the northern half of the system, and may in- 
dicate magnetic field shears or discontinuities, given that 
the non-thermal emission in galaxies is from synchrotron. 
The X-ray emission has a distinct minimum in the 'impact 
region', which is Hkely to be due to absorption of the soft 
photons (the blue area in fig. 2). Mid-infrared ISO data 
have shown that this is the regions of most intense star 
formation, including stars as massive as 60 Mq\ the ex- 
tinction in this region reaches Ay = 70 mag (Mirabel et 
al 1998). The blue area in fig. 2 can be seen trailing up 
to the northern nucleus (NGC 4038), following the region 
rich in dust lanes (see HST image in Whitmore et al 1999). 

All wavebands trace well the spiral arm that ap- 
pears to emerge on the SW side of the northern nucleus 
(NGC 4038), swings counter-clockwise and is easily trace- 
able for about 90 deg in radio and X-ray. At that point, 
NW of the nucleus, the radio arm dims but continues its 
swing for another 90 degrees, as does the Ha arm, and 
the HI emission. The diffuse X-ray emission instead is 
faint in this region. This region does not present any 
peculiarity in the mix of ages of stellar clusters (Zhang 
et al 2001), by comparison with other regions where dif- 
fuse X-ray emission is present, that would suggest a sig- 
nificantly older stellar population lacking the means of 
heating the ISM to X-ray temperatures. Hibbard et al 



(2001) point out that the HI kinematics in this area sug- 
gest infah of HI onto NGC 4038. It may be possible that 
the infalling cold gas interact with and displace the hot 
gas. A general pattern of avoidance between cold and 
hot ISM can be seen for example in the merger remnant 
NGC 5128 (Karovska et al 2002). From the velocity maps 
and Nh in Hibbard et al (2001), we estimate an infall ve- 
locity for the HI cloud ~ 100 km/s, and a mass density 
/9 ~ 1.6 X 10~^^ g cm^'^, corresponding to 1 H atom per cu- 
bic centimeter. The ram pressure of this cloud would then 
be p ~ 1.6 X 10~^° dyn cm~^, larger than the typical 
thermal pressure of the hot ISM (see Table 7) . Our lowest 
estimate of the ram pressure is '--^ 3 x lO^^^dyn cm^"*, com- 
parable to the pressure estimates for the hot ISM. How- 
ever, the pressures in Table 7 are derived for the most 
luminous X-ray regions, so are an overestimate of the pres- 
sure in the region of the northern spiral arm. We conclude 
that displacement of the hot ISM by infalling cold neutral 
hydrogen cloud is possible. 

Interaction between cold and hot ISM is also suggested 
by the relative distribution of CO and X-ray emission. In 
the impact region, X-ray colors and spectra may suggest 
extinction as the cause of the low X-ray surface bright- 
ness, although, given the uncertainties in the fits, a situa- 
tion with an embedded obscured hard component, and a 
'surface' low-emission measure soft component, is also pos- 
sible. However, the extinction interpretation would be in 
agreement with emission originating from a region of star 
formation deeply embedded in the molecular clouds, since 
the CO emission peaks in this area. Other CO/X-ray over- 
lap or avoidance areas are not immediately ascribable to 
extinction. These regions include the X-ray gap in the NW 
spiral arm. and regions where the X-ray emission is instead 
intense: the 'superbubble' in the NW arm discussed in Pa- 
per I (regions Rl and R2 of fig. 3b), and the two nuclear 
regions. As we have shown in §5.1, the pressures in the 
cold (CO) and hot (X-ray) ISM are similar and suggest a 
situation where cold and hot ISM may be commingled and 
in equilibrium. Near the most intense region of diffuse X- 
ray emission at the southern nucleus (NGC 4039) there is 
a marked minimum of both HI and CO. This may suggest 
a hot nuclear outflow pushing aside the cold ISM, reminis- 
cent of the morphology of M82 and NGC 253 (Watson et 
al 1984, Fabbiano & Trinchieri 1984). 

6. SUMMARY AND CONCLUSIONS 

This paper presents a first look at the complex hot ISM 
of the Antennae galaxies, based on the first Chandra ob- 
servation of this merging system. 

After producing an X-ray-color image, that gives us an 
immediate visual representation of the spectral complex- 
ity of the diffuse emission, we have obtained a quantitative 
estimate of this complexity by means of spectral fits. At 
least two temperatures (with different Nh) are needed to 
fit the most luminous regions, including the two nuclei, 
and what appear to be luminous superbubbles (Paper I) 
in the N-W star-forming arm of NGC 4038. The two tem- 
peratures are in the vicinities of 0.3 and 0.7 keV, with the 
hotter component being more highly absorbed. However, 
the fits (even adding a third component) are still not a 
good representation of the data. We have explored the 
effect of non-solar abundances in these fits. While the 
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fits improve for very low abundance values (as already re- 
ported with ASCA, Sanson! ct al 1996), they are still not 
good. Moreover, given the clear complexity of the emis- 
sion field, suggested both by the X-ray-color image and by 
the spectral results, spuriously low value of the abundance 
may result because of a simplistic choice of fit models (e.g. 
Weaver et al 2000; Kim & Fabbiano 2003). 

From these results we derive estimates of the hot ISM 
parameters for these three regions. We derive cooling 
times of ^ 10^ — 10^ yr and masses of hot ISM of ^ 
10^ — 10^ Mq for each region. In the two nuclei, these 
masses are much smaller than those present in molecular 
gas (Wilson et al 2000), but the thermal pressures of the 
two phases of the ISM arc comparable, suggesting equi- 
librium. However, this conclusion depends on crucial as- 
sumption on the spatial distribution and filling factor of 
the ISM. It is also possible that, at least for the hottest 
component of the hot ISM, the thermal pressure may ex- 
ceed that of the cold medium. For the two nuclear regions, 
where a comparison is possible, the supernova rates are of 
the order of magnitude of those estimated by Neff & Ul- 
vestad (2000) for the same regions on the basis of radio 
continuum observations at 6cm. 

We have compared the distribution of the diffuse X-ray 
emission with that of the 6cm radio continuum, Ha (us- 
ing an archival HST WFPC-2 image), HI (Hibbard et al 
2001) and CO (Wilson et al 2000). While there is a gen- 
eral resemblance between difii'erent wavebands, that may 



be related to a common link with star formation phenom- 
ena (Zhang et al 2001), this comparison shows significant 
difi^erences. In particular, it shows very clearly the effect 
of obscuration on the X-ray emission, that has a minimum 
(where only hard photons escape) coincident with the large 
radio continuum and CO peaks in the impact region of the 
two galaxies. We also note a displacement of the diffuse 
X-ray emission on the northen side of the system, rela- 
tive to the spiral arm visible in Ha, radio continuum and 
HI. The HI velocities suggest infall in this region (Hibbard 
et al 2001), that results in a ram pressure exceeding the 
thermal pressure of the hot gas. This could explain the 
displacement of the X-ray emission towards the south. 

These results demonstrate the richness of the ISM of the 
Antennae galaxies and justify further in dept-studies. A 
much deeper X-ray image is being collected with Chandra 
that will allow in time far more detailed studies of indi- 
vidual regions. 

We thank the CXC DS and SDS teams for their efforts 
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the reduction (SDP) and analysis (CIAO). We thank our 
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ital versions of their data that were used in this work, and 
D. Burke for help in producing fig. 6. The archival Chan- 
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was supported by NASA contract NAS 8-39073 (CXC). 
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Fig. 1. — Multi-color image of the Antennae prior to point source removal. Point sources to be removed are indicated by white ellipses. 
The color scale has been modified to enhance faint features. 



Fig. 2. — Multi-color image after point source removal. The color scale has been modified to enhance faint features. 



Fig. 3. — (a) Raw full band image of the Antennae after extracting the point sources indicated by their ellipses and numbered following the 
notation of Paper II. (b) Adaptively smoothed (0.3-6.0) keV band image of the diffuse emission of the Antennae, showing the regions used to 
extract the X-ray spectra. 



Fig. 4. — Confidence contours for the RS-I-(RS) fits for the regions shown if Fig. lb. Red contours corresponbd to the unabsorbed 
temperature component while blue contours correspond to the absorbed component. The crosses mark the best fit parameters. 




Fig. 5. — Spectra of the Northern and Southern nucleus and Regions 1 and 2 (see text for details) together with the best fit solar abundance 
RS-I-(RS) model and the residuals in a. The dashed and dot-dashed lines show the two thermal components. 
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Fig. 6. — Image of the low-filtered diffuse radio emission in NGC4038/4039 at 6cm(left) and at 20cm(riglit). Both images have been restored 
with a beam 6.59" x 4.25", and smoothed with a 7.6" x 7.6"low-pass filter to remove the compact radio sources. The 20cm image has an rms 
of 52mjy/beam, the 6cm image has an rms of 18mjy/beam. 



Fig. 7. — Gray-scale representation of the diffuse emission's 20-6cin spectral index distribution in NGC4038/9. This image shows only 
regions where both input images have a flux density of at least lOcr. The majority of the diffuse emission has a spectral index a <-0.5 (Si/ cx: v"") 
implying non-thermal synchrotron emission from relativistic electrons. Three regions have a. > —0.5; all are locations of significant current 
star-formation. Filamentary structures in the non-thermal emission probably indicates shocks and/or regions of magnetic field compression. 



Fig. 8. — Radio/X-ray image of the diffuse emission in the Antennae. The radio image (green) is from the 6 cm low-pass filtered image 
in fig. 6. The X-ray image (red) is from the (0.3-6.0) keV band and is adaptively smoothed. The countours give sketchy outlines of the Ha 
emission from the HST WFPC2 data (see Paper I), and identify the most active HII regions. The circles in the bottom right of the image 
indicate the resolution of the radio and X-ray bands. In the case of the X-ray band they indicate the minimum and the maximum smoothing 
scales. 



Fig. 9. — A combined HI (red), CO (green), and X-ray (blue) image of the Antennae. The HI data arc from the maps of Hibbard et 
al. (2001) and the CO data are from Wilson et al. (2001). The X-ray data are in the (0.3-6.0) keV band. The circles in the bottom right of 
the figure indicate the resolution of each image. For the X-ray band they indicate the minimum and the maximum smoothing scales. 
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Table 1 

Spectral Fits for the Temperature Map (Z = Zq) 



Region 
(1) 


Net counts 
(2) 


(3) 




(4) 


(5) 


(6) 


EM1/EM2 
(7) 


(dof) 
(8) 


1-A« 


33.4 ± 13.1 


0.22tl-Xt 


0.1(< 0.2) 








5.4 (8) 


l-B" 


135.1 ± 16.7 


0.601°:^, 





14+0.06 
^^-0.07 








13.3 (14) 


1-C 


285.4 ±20.6 


0.034(< 0.1) 





„o+0.03 
^*'-0.04 


< 0.17 


5.74t",:f, 


0.52 


22.6 (19) 


1-D 


247.4 ± 19.4 


0.034(< 0.09) 





on+004 

''^-0.04 


0-30l°j^ 


o.7oi«:;° 


0.79 


9.5 (17) 


1-E« 


130.2 ± 16.6 


^ ■ ' —0.19 





14+0.03 

-^^-0.08 








15.0 (14) 


2-A 


207.9 ± 18.5 


94+" ■''3 

— 0.26 





Qg + 0.03 

^^-0.02 


> 3.19 


6.6(> 0.7) 


12300 


11.1 (15) 


2-B 


767.0 ± 30.1 


0.034(< 0.063) 





OR+0.03 
"^"-0.04 




69+° °'' 


0.37 


44.3 (38) 


2-C 


1520.7 ±40.7 


0.0181°;°?^ 





qi +0.20 
•^^-0.02 






0.25 


123.0 (58) 


2-D 


554.4 ±26.4 


0.034(< 0.087) 





9(,+0.02 
^"-0.03 




0-85t°:l| 


0.60 


44.1 (31) 


2-E 


549.7 ±26.4 


0.034(< 0.056) 





9C+0.04 
^"-0.03 


0.37±S;f, 


Q y9 + 0.28 
^- '^^-0.07 


0.89 


48.7 (29) 


3-A 


928.2 ± 32.5 


0.034(< 0.054) 





9C+0.05 
''"-0.03 


0.4i;°;f, 


0.06 


0.39 


44.7 (46) 


3-B 


866.5 ± 31.4 







9r+0.03 
''"'-0.03 




^- ' ^ — 0.04 


0.15 


54.7 (42) 


3-C 


751.0 ± 29.9 


0.79l°:fg 





,,+0.03 
^■^ — 0.02 


> 0.12 




209 


39.9 (41) 


3-D 


694.6 ± 28.6 


0.034(< 0.049) 





9r+0.03 
^^-■-0.03 


0.44tS;;| 


0.69t°r6 


0.35 


40.6 (36) 


3-E 


794.7 ±30.5 


0.034(< 0.052) 





9, +0.03 
'^-'--0.02 


o.24iS:24 


0.69l°:J? 


0.77 


58.9 (42) 


4-A« 


202.0 ± 18.3 


0.8llo'22 





ir+0.14 
^'-*-0.22 








15.1 (17) 


4-B 


368.0 ± 22.1 


0.034(< 0.11) 





on+0.05 

•^^-0.06 


o.8otg:f. 


0.88t°„:l? 


0.09 


23.2 (23) 


4-C 


483.1 ± 24.9 


0.06+°°S 

— 0.02 





9fl+0.04 
'^"-0.04 


0.26+°-?! 

— 0.2b 


0.791°°;^ 


0.58 


16.6 (26) 


4-D 


264.6 ± 20.3 


0.034(< 0.94) 





9fi+0.05 
^°-0.04 


"■^-^-0.41 


(, 70+0.22 


0.60 


12.6 (18) 


4-E 


127.2 ± 16.5 


0.034(< 0.1) 





90 + 0.03 
^°-0.03 


99+''^-" 


> 2.24 


0.33 


11.3 (11) 


5-A 


141.7 ± 17.0 


1 4+1-1 





(,7+0.05 
-0.03 


a .<t: + 4.45 


0.351^0°,? 


654.4 


8.3 (11) 


5-B 


714.4 ±29.1 


0.034(< 0.079) 





of,+0.03 


f) ,7+0.13 

"•^'-0.18 


0.83l°:°o| 


0.28 


46.5 (38) 


5-C 


2101.0 ±47.3 


0.034(< 0.047) 





„n+0.02 
•^^-0.02 


42+° i° 
"•*-^-o.io 


n 70+0.05 

'^•'8-0.05 


0.32 


104.5 (68) 


5-D 


374.2 ± 22.7 


0.034(< 0.06) 





,n+o.05 
■J^-0.12 


< 0.29 


' -^-0.13 


1.76 


28.0 (22) 


5-E 


131.9 ± 16.5 







90+0.12 

'^"-0.10 


> 17.9 




<0.0001 


11.2 (11) 


6-A« 


116.2 ± 16.1 


0.034(< 0.126) 





9-.+0.03 
^ ' -0.05 








11.7 (13) 


6-B' 


150.6 ± 17.1 


n QQ+O-22 





19+0.03 
^^^-0.03 








14.1 (15) 


6-C 


293.9 ± 20.7 


U'-^«5_0.20 





, 0+0. 06 
^^-0.06 


< 0.29 


0.48l°o:f7 


3.90 


15.5 (18) 


6-D 


139.7 ± 16.8 


0.034(< 0.15) 





,0+0.11 
•^^-0.07 


< 0.52 


l-26l°-7,^ 


0.52 


8.2 (11) 



'^'^ Absorbing column density for total spectrum (a) and the second thermal component (c) in units of 
10^^ cm^^. Quoted errors are at the 90% confidence level for one interesting parameter. 

^''^ Temperature for the unabsorbcd (b) and the absorbed (d) thermal component in keV. Quoted errors are at 
the 90% confidence level for one interesting parameter. 

^ Fits with only one thermal component. 
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Table 2 

Spectral Fits for the Temperature Map (free Abundance) 



Region 

(1) 


(2) 


kTj 

(3) 




Z" 

(4) 


(5) 


(6) 


(7) 


IP A/T / TP Tl/r 

Cj lVl\ / Hi lVi2 

(8) 


1,2 

X (dot) 
(9) 


1-AS 




U.UO_(j 03 


0.02(< 5.0) 










3.3/7 


l-B" 


0.24(< 0.678) 


22+°-^° 


U.UO_o 021 










9.7/13 


1-C 




r, OO + 0.05 
U.^»_0,08 


0, 


„,+4.96 
■"*-0.03 


< 9.02 


64.0'' 


5.00'' 


19.02 


20.0/17 


1-D 




21 + °° = 
"■^^-0.05 


0, 


24+4.76 
■^^-0.19 


n , (5+0. 30 
0-3610.36 


Q 70+0.15 

"• '" — 0.13 


4.00(> 0.26) 


6.96 


9.1/15 


1-E9 


0.19(< 1.05) 


0.03(< 0.24) 


0, 


90 + 0-23 
■^^-0-19 










11.2/12 


2-A 


74+0 39 

'^-0.45 


ll+° i3 

Uii_0.04 


0, 


nQ + 0-42 

■"''-o.os 


'J-oa-0.85 


4.48(> 0.45) 


5.00'' 


9012 


7.8/13 


2-B 




n 97+0.06 


0, 


,7+4.73 
■^'-0.18 


U..3U_o 26 


70+0.07 

"• '" — 0.07 


0.45(> 0.14) 


0.70 


42.5/36 


2-C 


f) (,0+0.05 


n O1+0.03 
"■•'^-0.03 


0, 


o,+4.67 
■'^''-0.22 


(, 9C+0.31 
"-^0-0.14 


o-7it°:°oI 


09+°-^* 

"-""-0.04 


0.24 


111.5/56 


2-D 




n ,7+0.03 
"■^'-0.05 


4, 


qO+0.10 
■"^-4.52 


02+°-^" 
"-"^-0.02 


0.89t°o:l° 


i4+°-20 

"-^*-0.09 


0.05 


38.9/29 


2-E 




"■■^"-0.03 


0, 


9^ +0.41 
■^-^-0.13 


1 90+^°^ 


09+° °* 

"•""-0.01 


4 KK + 0.45 


0.0001 


37.3/27 


3-A 


06+°-°" 

U.UD_Q Q2 


27+° °'^ 

U-^ ' -0.04 


0, 


1 , + 0.25 
■^•^-0.07 


40+°-^^ 

"■^"-0.34 


69+°-^° 

"-""-0.06 


1 29+^-^^ 

^-^"-1.10 


3.54 


38.3/44 


3-B 


n i3+° i5 

^•^•5 — 0.09 


r. or+o.oe 

"•■'U-0.06 


0, 


(,7+0.41 
■"'-0.05 


64+°-^" 

"-"^—0.27 


72+0.07 

"• ''^-0.05 


1 1 9+3.88 
^•^''-0.86 


3.00 


51.1/40 


3-C 


29+°-^° 


22+° °'' 


0, 


r,9+4.98 
■"^-0.02 


< 0.36 


77+°-" 

"- ' ' -0.10 


1 fj+1.38 
"• J-O — 0.10 


14.68 


25.8/39 


3-D 


f) „,+0.04 


0.30i°o°« 


0, 


.. +0.10 
■^^-0.05 


"-^''-0.32 


70+°-°s> 

"- '"-0.06 


4.99(> 0.39) 


10.52 


28.8/34 


3-E 


f, po + 0.06 


0.23t°°ot 


1, 


,o + 3.62 
■■^"-1-28 


< 0.35 


71+0.09 

"■ ' ^ — 0.08 


13+1-49 
"■l''-0.05 


0.11 


45.1/40 


4-Af 


n ^^7+1-13 

^■•^'-0.19 


27+° " 

"■^'-0.20 


0, 


05+°-" 
■"°-0.04 










11.6/16 


4-B 




,c+0.24 
"••5u_o.09 


0, 


(,0+0.69 
■""-0.07 


37+0.37 
"-"' -0.55 


91+°-^i 

"•"^-0.14 


45+4-55 

"•^^ — 0.36 


0.52 


19.7/21 


4-C 




0.26l°o°o? 


0, 


(,7+4.93 
■"'-0.05 


20+0-37 

"•^"-0.20 


79+0-14 

"• '"-0.11 


1 42+3-58 

^•^■'-1.31 


16.47 


14.0/24 


4-D 




0.29i°:°^7 


0, 


9n+4.80 
■^^-0.15 


52+°-=^ 


0.79lS:?5 


4 go+o.oi 

''-»»_4-74 


10.34 


11.1/16 


4-E 




■?n+°-^° 


0, 


1 , + 4.87 
■-'-■'-0-10 




1.44(> 0.64) 


5.00'' 


1.09 


10.0/9 


5-A 


1 12+^-" 


0.08lS:o4 


0, 


p, + 4.97 
■"'^-0.03 


1 Q+33.8 
J^"-y-10.9 


"-^ ' -0.20 


03+''-°'' 
U.UO_o 03 


5.10 


6.2/9 


5-B 




,0+0.08 
"••5 -^-0.07 


0, 


.., + 1.18 
■^^-0.09 


40+°-^^ 

"•^" — 0.40 


0.83jlo;j9 


, ,,5+1.64 
O.OD_3 07 


7.21 


44.1/36 


5-C 


„ f,,+0.06 


34+0 04 

^•^^-0.04 


0, 


■^^-0.30 


13+0-14 

"• J^^-0.13 


n S5O+0.06 
"•°''-0.06 


f, ,Q+0.19 
0.i»_o.o6 


0.25 


84.8/66 


5-D 


05+°-°** 


42+° i° 


0, 


1 , + 0.11 

■-''^-0-05 


89 1+*^-° 

'''--'-83-6 


i7+°-°4 

-0-04 


4.83'' 


<0.0001 


21.9/20 


5-E 




"•'■^8-0.10 


0, 


,g + 4-85 

■-'°-0-12 


-1 ,o +27.0 
ll = --45.7 


,7+0.08 

"- ' —0.03 


4.32'' 


<0.0001 


10.5/9 


e-A" 


0.034(< 0.16) 


■?i+o io 

'J.ol_o 07 


0, 


gg + 0-18 

■"°-0.04 










6.4/12 


e-B" 


n 7Q+0.38 


f, jo+O lO 
"•^•^-0.03 




0.04'' 










6.4/12 


6-C 


14+°" 
^•^^—0.09 


o.4ii°o:il 


0, 


■in+0.09 

■■'"-0.04 


10 8+''-^° 


63.9'' 


0.88'' 


5.30 


18.1/16 


6-D 


0.03(0.36) 


0.36t°o:f3 




0.12'' 


< 1.18 


,1+5.98 
"••^1-0.75 


5.00'' 


15.79 


7.2/9 



^'"^ Absorbing column density for total spectrum (a) and the second thermal coinponcnt (c) in units of 10^^ cm ^ . Quoted 
errors are at the 90% confidence level for one interesting parameter. 

^'■^ Temperature for the unabsorbed (b) and the absorbed (d) thermal component in kcV. Quoted errors are at the 90% confidence 
level for one interesting parameter. 

Abundance for the unabsorbed (b) and the absorbed (d) thermal component in keV. Quoted errors are at the 90% confidence 
level for one interesting parameter. 

^ Fits with only one thermal component. 

^ Unconstrained parameter (lo" range greater than allowed range for the parameter). 
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Table 3 

Spectral Fits for Selected Regions (Z=Z©) 



Region 
Net counts ± error 



(1) 



RA (J2000) 
Dec (J2000) 
Area (sq. arcsec) 

(2) 



RS+(PO) 
Nn.tot" kT? 
Nh^ r 

(dof) 
(3) (4) 



RS+(RS) 
NH.tot" kT? 

(dof) 
(5) (6) 



North Nucleus 


12: 


01: 52.8 


0.034(< 0.2) 


n cc+0.06 
U.DO_o.o8 


0.034(< 0.045) 


n ct;+0.03 

U.DO_Q Qg 


604.7 ±25.4 


-18: 


52: 08.5 
222 


0.15(< 0.28) 


^•^-0.5 

38.9 (29) 


< 0.015 


52.3 (29) 


South Nucleus 


12: 


01: 53.9 


0.034(< 0.058) 


n 7fi+o.03 

U.^D_Q 03 


0.034(< 0.045) 


n q9+0.02 
'J-<J^-0.03 


2180.4 ±47.5 


-18: 


53: 08.5 
488 


r. 99+0.03 
'-'•^'^-0.09 


q 4 + 0.7 

138.7 (67) 


41 


0.79tro! 
125.1 (67) 


Region 1 


12: 


01: 50.4 


0.034(< 0.08) 


fiti+o os 

U.00_o.o5 


0.034(< 0.05) 


20+0 0^ 


780.3 ± 28.9 


-18: 


52: 19.1 
235 


'-'•^^-0.09 


o qo+0.73 

64.9 (34) 


< 0.18 


n 70+O.O6 

58.9 (34) 


Region 2 


12: 


01: 50.6 


0.11(< 0.27) 


n qq+0.09 


0.034(< 0.08) 


97+0.04 


302.0 ± 18.0 


-18: 


52: 04.3 
137 


< 0.45 


^•°-0.8 

14.9 (13) 


"-"•^^^-0.36 


n 7^+0.30 

11.2 (13) 



'^•'^ Absorbing cohimn density for the total spectrum in units of 10^^ cm ^. Quoted errors are at the 

90% confidence level for one interesting parameter. 

^''^ Temperature of the thermal component in units of keV. Quoted errors are at the 90% confidence 
level for one interesting parameter. 
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Table 4 

Spectral Fits for Selected Regions (free Abundance) 



Region 



(1) 



Nn.tot 
(2) 



RS+(PO) 



r 



(dof) 
(3) 



RS+(RS) 

NH,tot" kT5 

(dof) 
(4) (5) 



Zf 
(6) 



North Nucleus 


0.14(< 0.23) 


0.63tr2l 


0.06(< 0.15) 


n O+0.08 

"••^-0.06 


17+1-89 




3.4(< 30.3) 


o 1 9+6.8 




0.71°^ 


12+019 
"•-i^-o.ig 




10+°°2 

"•■'-^-0.02 


37.5 (28) 


46+° ''^ 

^•^"^-0.31 


32.5 (27) 




South Nucleus 


in+o o2 

U.iU_o.o2 


74+0.03 

U-i4_0.07 


0.034(< 0.1) 


0.35t°;0« 


0.22(> 0.08) 




2.5(< 180.3) 


-0.45(< 6.5) 


17+013 


0.82tn^ 


r. 00+0.27 

'-'•^''-0.13 




n 14+005 


135.1 (66) 




99.3 (65) 




Region 1 


0.034(< 0.108) 


S7+°-" 

-0.07 


0.034(< 0.105) 




0.09(> 0.02) 






9.34(> 6.37) 


0.05(< 0.36) 


A 70+0.09 


0.29(> 0.09) 




o.()7^;i«,3 


41.8 (.'i'^) 




41.7 (32) 





'^''^ Absorbing column density for the total spectrum in units of 10^^ cm ^. Quoted errors 
are at the 90% confidence level for one interesting parameter. 

^'"^ Temperature of the thermal component in units of keV. Quoted errors are at the 90% 
confidence level for one interesting parameter. 



Table 5 

Spectral Fits for Selected Regions (Z=Z©) 



Region RSi + (RS2+PO) 

r (dof) 
(1) (2) (3) 



RSi + R.S2 

NH.tot" 

r 

(4) 



+ (PO) 

(dof) 
(5) 



RSi + (RS 

NH.tot' 

(6? 



2+RS3) 

kxf 
(dof) 
(7) 



RSi + RS2 

NH.tot"- 

(sf 



+ (RS3) 
kTj 
kT° 
(dof) 
(9) 



North Nucleus 


0.034(< 0.18) 


0.28lS:S^, 


0.034(< 0.21) 


0.28l°:J° 


0.034(< 0.72) 




0.034(< 0.075) 


f, OQ+0 05 




0.24(< 0.69) 


o.7ot°:ii 


0.27(< 0.83) 


0.72i°;i3 


f) 70+0.19 


0.651°-? 


68+° ™ 


8.1(> 2.7) 




9 Q7+1.08 


30.0 (27) 


2 96+^ " 


30.6 (27) 


16.2(> 1.9) 


33.2 (27) 


66+°-°° 


30.7 (27) 


South Nucleus 


0.034(< 0.10) 


„ oc+0.06 
U.OO_Q Q3 


0.034(< 0.15) 


n 07+0. 06 

"J.Oi _o.03 


0.034(< 0.072) 


-JK + O.OS 


0.053(< 0.073) 


0-36tS:S^3 




o.i5iS:i^ 


n 82+°-°^ 

>J.O.'_0.04 


0.33(< 0.6) 


0.85tS:S? 


0.03(< 0.08) 


n 43r+0.05 
0-85Io.o4 


< 0.09 


o.85i°o:;i 




2.9tgJ 


99.7 (65) 


-1.09 


101.2 (65) 


q+20 


112.8 (65) 




113.1 (65) 


Region 1 


0.034(< 0.13) 


^■^°-0.06 


0.06(< 0.14) 


n 2f; + 0.08 




o.o7i°:Si 


23+° l° 


o-06j:°g°,^2 




0.05(< 0.19) 


n 74+0 08 


< 0.32 




0.05(< 0.3) 


o.28iS:St 


0.032(< 0.29) 


0.281°;- 




3.53^Q 9 


42.4 (32) 


-0.87 


42.6 (32) 


' '^ — 0.08 


42.8 (32) 


79+"- 1° 


42.9 (32) 



^■''^ Absorbing column density for the total spectrum (a) and the hard component (c) in units of 10 cm . Quoted errors are at the 90% 
confidence level for one interesting parameter. 

^'"^ Temperature of the unabsorbed (b) and absorbed (d) thermal component in units of keV. Quoted errors are at the 90% confidence level 
for one interesting parameter. 
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Table 6 

Intensity of Spectral Components (Z=Z0) 



Region 


RSi + (RS2+PO) 


RSi -\~ RS2 


+(PO) 


RSi + (RS2+RS3) 


RSi + RS2 


+ (RS3) 




EM? (Lx) 


EM° (Lx) 


EM^ (Lx) 


EM^ (Lx) 


EM^" (Lx) 


EM^ (Lx) 


EM? (Lx) 


EM^ (Lx) 




Norm.'' (Lx) 




Norm.*^ (Lx) 




EM^ (Lx) 




EM^ (Lx) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


North Nucleus 


2.9 (0.38) 


5.3 (0.89) 


2.3 (0.30) 


2.6 (0.43) 


4.62 (0.60) 


28.8 (4.96) 


4.03 (0.52) 


4.87 (0.47) 




7.9 (4.53) 




9.6 (6.28) 




4.21 (0.40) 




22.8 (3.93) 




South Nucleus 


9.7 (1.32) 


26.1 (3.70) 


8.9 (1.24) 


17.1 (2.34) 


10.1 (1.39) 


21.2 (2.88) 


10.9 (1.52) 


20.3 (2.74) 




11.9 (7.38) 




23.7 (33.0) 




16.5 (1.27) 




16.6 (1.28) 




Region 1 


4.6 (0.60) 


5.8 (0.92) 


5.6 (0.73) 


5.4 (0.86) 


350.2 (338.7) 


17.4 (2.25) 


95.6 (96.6) 


19.7 (2.54) 




2.6 (3.95) 




2.4 (3.04) 




8.4 (1.26) 




7.95 (1.19) 




Region 2'' 


3.4 (0.44) 


6.2 (0.90) 






4.0 (0.53) 









2.4 (1.27) 



^'^ Emission Measure for the low (a) and high (b) temperature component in units of 10^ X (AttD^) cm ^ where D is the distance. In 
parentheses, absorption-corrected luminosity in the 0.1-lO.OkeV band, in units of 10^^ erg s~^ 

Normalization of the PO component (at 1 keV) in units of 10~^ photons/keV/cm^/s. In parentheses, absorption-corrected luminosity 
in the 0.1-lO.OkeV band, in units of 10^^ erg s~^ 

^ Emission Measure (and luminosity) for the third thermal component. The units arc the same as in the first two components. 

^ The parameters for Region 2 correspond to the RS+(PO) and RS+(RS) models instead of the RS+(RS+PO) and RS+(RS-i-RS) 
models respectively. 



Table 7 
Estimated ISM Parameters 



Region 


n 




Tc 


MjSM 


P 


SNR 




(10-2 cm-3) 


(10^^ ergs) 


(10^ yrs) 


(lO^Mo) 


IQ-^Myn cni-2 


(10-3 yrs-i) 


North Nucleus 


3.3i, 4.52 


0.5i, 1.62 


5i, 62 


0.3i, 0.42 


3.I1, 10.72 


0.2 




(1.5i, 2.02)« 


(1.2i, 3.62) 


(lOi, 132) 


(0.7i, 0.92) 


(1.4i, 4.82) 


(0.3) 


South Nucleus 


1.9i, 3.22 


3.O1, 142 


7.7i, 122 


1.7i, 32 


2.3i, 8.92 


1.2 




(0.9i, 1.42) 


(7.I1, 3O2) 


(I81, 272) 


(4i, 6.42) 


(l.li, 3.92) 


(3.0) 


Region 1 


2.2i, 2.52 


l.Oi, 3.52 


5.4i, 132 


0.7i, 0.82 


1.9i, 6.32 


0.3 




(l.Oi, I.I2) 


(2.I1, 7.72) 


(12i, 282) 


(I.61, 1.72) 


(O.81, 2.82) 


(0.7) 



^ Lower temperature RS component 
2 Higher temperature RS component 

^For an emitting cylinder of height 200 pc; values in parentheses are for a 1000 pc height cylinder 
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